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CN . Abstract 

. Dynamics of two anharmonic oscillators with interaction of the fourth order has been inves- 

["t I ! tigated. The conditions at realization of which system is integrable are established. The exact 

I analytical solution of the nonlinear equations in the case of adiabatic isolation of a system of 

PsJ ■ oscillators has been obtained. 

PACS: 76.60.-k 

^ ! In optics of weak light streams the polarization of an electromagnetic wave in substance does not 

■ depend on intensity of light and univocally corresponds to the polarization of emanation, impinging 
on the boundary of medium - vacuum. The situation cardinally varies in the case of nonlinear 
optics, when the factors of a refraction and absorption of substance become functions of intensity 
of a radiation [1]. At a particular combination of parameters of a radiation and nonlinear medium, 
the stationary polarization of light " on an inlet " becomes unstable " on an exit ", or there can 

, be a random modification of polarization in time [2-7]. The specified phenomenon is called as a 
O I spontaneous violation of polarization symmetry. 

^ ■ In a strong light field, when the nonlinear component of the electrical induction becomes notice- 

■ able with respect to the linear component, traditional exposition based on a material equation of 
^ • a medium, which represents a series with respect of terms of electric field is impossible. In works 

[8-9] the prospects of observation of amplitude instability for mediums in which nonlinear response 
is featured within the framework of a model of the Duffing oscillator were considered. In work [10] 
the nonlinear susceptibility of the ensemble of chaotically oriented mirror - asymmetrical molecules 
has been estimated. Assumption that each of molecules can be represented by the nonlinear system 
of interacting oscillators in the field of electromagnetic wave has been done: 

X + uj'ix + ay + Tx + ax^ + b(y'^ + 3x^y) + cxy"^ = — - expiiKz — ) 

m 2 

X + ujly + ax + Vy + ay^ + b{x^ + Sy'^x) + cyx"^ = — - exp^iK^ — ) (1) 

m 2 

Where ujq is a natural frequency of small oscillations of oscillators creating a molecule, e and m 
electrical charge and mass accordingly, Ez, Ey projections of amplitude of an electric field, D is 
a distance between oscillators, a , b , c are the anharmonic constants naturally of oscillators and 
connection between them, T is the parameter describing attenuation. With the help of summation 
of the dipole moment of molecules 

dx,y = ^ ( ^ ) exp{±iKz^) 

and averaging with respect of small volume, in the work [10], the expression for macroscopic 
polarization was obtained, and was shown that in mediums with nonlinear gyrotropy the polarization 
instability takes place. 
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Despite of it, due to the great importance for the nonhnear optics, the immediate research of the 
system of nonhnear oscillators (1) is of interest. 

The purpose of given work is the qualitative, theoretical research of the equations (1) and estab- 
lishment of relations between parameters of a system at which equation (1) will be integrable. 

Let us consider at first the case of weak nonlinearity. We also suppose that the amplitude of the 
external field is small: 

^0 » 

m 

Considering terms with cubic nonlinearity as a perturbation, let us rewrite system (1) in following 
form: 

X + uj^x + ay^ -l^f{x, x, y, y), 

y + Uoy + ax^ -lig(x, x,y,y), (2) 



where 



f{x, X, y, y)^x^ + y^ + Sx'^y + xy'^ + -x 

r 

9{x, X, y, y)^x^ + y^ + Sxy'^ + x^y + -y (3) 

A* 



here are used the following designations: 
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and transformation to the dimensionless variables and parameters (where - a^ means a Bohr radius) 
has been done. For /j, — the set of equations (2) has the solution: 

x = a-^ sin(A;^t + + sm{k^t + jS^) 
y = —a^ sin{k^t + 13^) + si'ii{kj, + /5J (4) 



Here q;^, cKj, (5r^ are the constants of integration, — \JuJo — a^, k^ — ^Jujq — 

Solution of the set of equations (2) for 7^ we shall search in the form (4), considering that 
ttj , 0:2 , A ) 1^2 s-re slowly varying functions of time. The additional conditions which we shall impose 
on the functions ,a^,j3^ , {3^ for their determinancy consist in the following: first order temporary 
derivatives of x and y are having the same form, as in the case of constant , ctj , /J^ , /Jj- As a result 
we get: 

sin{k^t + /3J + sin(A;2t + 13^) + a^jS^ cos{k^t + /3J + a.^(3.^ cos{k.^t + (3.^ = 
— dj sin(A;^t + /3J + dj sin(A;2t + P^) ~ <^iA cos{k^t + /SJ + 0-2(^2 cos(A;2t + (d^) = (5) 

Calculating temporary derivatives of x, and y and substituting obtained expressions in the 
equations (2) one can get: 

djfcj cos(A;^t + /9J + d2A;2 cos^k^t + P^) — k^a^P^ sm{kj. + /SJ + k^a^(3^ sm.{k^t + 13^) = fif* 
—a^k^cos{k^t + /3J + a^k^ cos{k^t + P^) — k^a^jS^ sin{k^t + /3J — k^a^(3^ sm{k^t + 13^) = fig* (6) 



Where 



/* = /[«! sin(A;ii + + 0:2 sin(A;2* + ^2) ; cos{kit + + a2k2 cos{k2t + (32) ; 

-ai sm{kit + + 0:2 sm{k2t + P2) ; -aiki cos{kit + + q;2A;2 cos(A;2^ + P2) ], 
g* = g[ai sm{kit + Pi) + 0:2 sm{k2t + P2) ; oiiki cos{kit + Pi) + a2k2 cos(/c2^ + P2) ; 
—ai sm{kit + Pi) + a2 sin(A;2^ + P2) ; —o-iki cos{kit + Pi) + a2k2 cos{k2t + P2) ] 

The equations (5), (6), (7) are a set of equations for the determination of q;^ , ctj , /3j , P^ 
After solving of these equations with respect to d^one can obtain: 



dai 
~df ~ 

sin ^ sin 77 

— sin ^ sin rj 

ki cos ^ /C2 cos 1] 

-/cicos^ k2 con rj 



Ai 
A 

cos ^ COS 77 

— COS ^ COS t] 

—kisin^ —k2sinr] 

A'Jisin^ — A;2sinr7 



sin T] COS ^ 

sin Tj — cos ^ 

A*/* ^2 cos 77 — /ci sin ^ 

/ig* k2 cos 77 /ci sin ^ 



COST] 

cost; 
-k2sinr] 
-k2 sin 77 



(7) 



here ^ = /cit + /3iand 1] = k2t + /32. 

Making similar operations with respect to ai ,a2 , Pi , P2 after unwieldy calculations we get: 
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da2 

~dr ^ k2{kl - kl) 
dpi 11 
dt ^ ki{kl-kl) 

dP2 _ fjL 

It ~ ~ki{kl - kl) 



[af* + ag*] cos 77, 
[af - ag*\ sin^, 
[af* + ag*] sin 77. 



(8) 



The obtained equations (8) represent a set of equations (2), transformed to another variables. 
Let's assume, variation of the cki , 0:2 , /3i , P2 are slow than oscillations in the initial dynamical system. 
Averaging the obtained equations with respect to phase |^ and |r for ai , q;2 , /3i , P2 we shall get: 
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Pi 



2a^ 




— a 








2a ^ 
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(Fi-Gi), 
{F2 + G2), 
(Fs-G,), 
iF, + G,). 



(9) 
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Here 

^1 = ^ / / /* ^1 = ^ / / ^* ^^^^^ 

ivr^ Jo JO /vr"^ Jo Jo 

]^ /■27r i'2tt \ r2TT /•27r 

^^^2n^Jo Jo -f*^^^^^^^^''^^^ Jq Jq 9* cos rjd^dr) 

/ / f* si'^^d^dr];G3= — / g* sin ^d^dr] 

Jo Jo ZTT^ Jo Jo 



(10) 



]^ /'27r /'27r ]^ r27T r27r 

■^'^^27r2Jo io ^*^^^^^^^^''^^^ 2tt^ Jq 7o 9* sin rjd^drj; 



The right sides of the equations (8) do not depend on (3i and /?2- Therefore they serve for the 
investigation of ai and a2- Substituting (7) in (10) and integrating finally we get: 
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ai{t)^ai{0)e ^H^^ , ai{t) ^ ai{0)e ^(^F^ . (11) 

It is easy to see from (11) that in case of weak nonlinearity, the nonlinear terms does not affect 
on the form of motion. Therefore the application of the method of slowly varying amphtudes in 
case of strong nonhnearity made in [9] is incorrect. As a result in the system originates damping 
oscillations and the phase diagram represents a steady focal point. This outcome can be considered 
as a singularity of the system (1). These reasonings are proved also by numerical calculations (see 
Fig. 1). 

Fig. 1 Phase diagram obtained by numerical calculations for the values of the dimensionless 
parameters: Uq = 10 , a = 1, F = 0, 5 , / = 0, 1 , a = 0,6, 6 = 0,5, c = 0, 5 

As show numerical calculations at the action on the system of a strong exterior field 

2 ^ 0.1^ (12) 

mao 

in the phase diagram there is a limit cycle [11-13] (sec Fig. 2). 

Fig. 2 Phase diagram obtained by numerical calculations for the values of the dimensionless 
parameters: 

u;2 ^ 10, a = 1, F = 0.5, / = 1.7, a = 0.6 6 = 0,5, c = 0.5. 
In case of strong nonlinearity the application of the method of slowly varying amplitudes becomes 
impossible. However below will be shown that under particular conditions is possible to obtain 
analytical solutions of the set of equations (1) even in the presence of interaction terms of the fourth 
order. 

It is easy to see that in case of absence external field and damping, (i.e. in case of adiabatic 
apartness) the equations (1) can be generated by the Hamiltonian 

H = + + 00^/2 (x' + y") + axy + a/4{x^ + y^) + c/2 x'y' + 6(xy=' + x^y) (13) 

After the transformation to the variables qi = x + y,q2 = x — y and supposing that c = 3a one 
can obtain: 

= f + + + (14) 

H2 = § + ^^4 + -yo-a)ql 

The similar result can be obtained with the using of more general method (the method offered 
by P. Lax [14]). This method allows solving the problem even in case when determinations of new 
variables are not so trivial. 



Let's initiate the analysis of the Hamiltonian (13). The corresponding equations look like 



Ql,2 

1 d , 



a + 6 3 uJq ±a 



2 ^1'^ 2 

Representing q^ ^ as g\ j = 2^^(^i,2)^^^'i substituting in (15) we get: 

^/\a + b 



(15) 



?1,2 



(16) 



where 



0,. - 4,A0) + ^(?,.(o))^ + ^ (?,,(o))^ 

o' ± a + ^J{u;^±ay-ie^,2\a + b\ 
\a + b\ ' 



± a - ^^ {uj^ ± ay - 4^1,2 |a + 6 
\a + b\ 



In (16) we have made assumption that a + b < because integration of (17) can be done only in 
this case. The integration of (16) gives 



1 p. . QlAt) C^ 

F(arcsm — ; -^-ri-) 



a + b\t 



$2 

1.2 



(17) 



where F(...) is the elliptic integral of the first kind [15]. Having produced a reversion of the 
expression (17) finally we get: 



?l,2(^) = 0i,2S^( 



^4 t-<' 

O ^1,2'' ) 



$2 

1,2 



(18) 



where sn (....) is the elliptic sine of the Jacobi [15]. 

The expressions (18) are general solutions of the equations corresponding to the Hamiltonian 
(13). By means of inverse transformations x — ^li^ ; y — they allow to express solution of the 

initial problem in the analytical form through higher transcendental functions. In case of the initial 
conditions , = \ ^^^^ 

the solutions become considerably simpler [15] and takes form of solitary waves [16] 
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'^° + "tanh,r^t + 



\a + b\ 



^°-"tanhJ^t) 



\ \a + b\ 



^°' + "tanhJ^t- 



\a + b\ 



tanh 1 ^^t) 



\ \a + b\ 



(19) 



For obtaining of solutions (20), (21) wc neglected a dissipation of energy and influence of external 
field, i.e. we suppose that the system (1) is adiabatic isolated and change of energy during the one 
period of oscillations can be neglected. Taking into account that the period of the solution (18) is 
defined by the expression 



1,2 

(where K (...) is the complete elhptic integral [15]), it is possible to announce that expressions 
(18) and (19) are valid for the intervals of time 



1,2 

In case of intervals of time t > "^''^"^" -T, under the action on the system of external constant field 
the obtaining of the analytical solution is impossible and it is necessary to use numerical methods. 

In case of strong nonlinearity for the particular values of parameters in the system can occur 
chaotic oscillations [11-13] (see Fig 3, Fig 4, Fig 5). 

Fig. 3 x{t) as a function of time obtained by numerical calculations for the values of the dimen- 
sionless parameters: o'q = 10, a = 1, T = 0.1, / = 1.7, a = 3, b = 2.5 , c = 2.5 

Fig. 4 as a function of time obtained by numerical calculations for the values of the dimensionless 
parameters: a;Q = 10, a = 1.2, F = 0.1, / = 1.5 , a = 3, 6 = 2.5, c = 2.5 

Fig. 5 Projection of trajectory of phase point motion on the space x{t), dx{t)/dt, y{t) for the 
values of parameters: a;Q = 10, a — 1, V — 0.1, f — 1.7 , a — 3, b — 2.5 , c — 2.5 obtained 
by numerical calculations. 

It is ease to see from Fig. 3, Fig. 4 and Fig. 5 that motion is realized in the restricted region of 
phase space. The external field compensates energy dissipation and we have not focal points. On the 
other hand motion is not periodically and as a result trajectory is not closed (see Fig. 5). In another 
words it means that ratios between frequencies of oscillations of x{t), x{t), y{t) are not rational 
numbers. So in case of strong nonlinearity, under the action of external constant field in the system 
occurs chaotic non-periodical oscillations and system is not integrable. 



(20) 
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